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ABSTRACT 

The similarity between observed velocity structures of Al III and singly ionized species 
in damped Lya systems (DLAs) suggests the presence of ionized gas in the regions where 
most metal absorption lines are formed. A simplified model consisting of Region 1) a 
plane-parallel ionization bounded region illuminated by an internal radiation field, and 
Region 2) a neutral region with a negligible metal content is considered. We calculate 
photoionization equilibrium models for region 1, and constrain the ionization parameter 
by the observed iV(Al III)/A r (Si II) column density ratio. Under these conditions we 
find that ionization effects are important. 

If these effects are taken into account, the element abundance ratios in DLAs are 
quite consistent with those observed in Milky Way stars and in H II regions of local 
low-metallicity blue compact dwarf galaxies. In particular we cannot exclude the same 
primary N origin in both DLAs and metal-poor galaxies. No depletion of heavy elements 
on dust grains needs to be invoked, although our models do not exclude the presence of 
little depletion. 

Although highly simplified and relying on the strong assumption of a significantly 
lower metal content in region 2, our model appears to be supported by recent data on 
a local DLA and it is not in contradiction with the current knowledge on high redshift 
DLAs. If correct, it offers a clear simplification in the understanding of heavy element 
abundance ratios in DLAs and their comparison with the local Universe. 

Subject headings: galaxies: abundances — intergalactic medium — quasars: absorption 
lines 
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1. INTRODUCTION 

Damped Lyman a systems (DLAs) are the class of 
QSO absorption systems showing the largest column 
density, with 7V(H I) > 2xl0 20 cm" 2 . They are con- 
sidered to be the high-redshift progenitors of present- 
day luminous galaxies (e.g. Wolfe 1990) and dominate 
the neutral gas content of the Universe (Wolfe et al. 
1995), a mass comparable with that of the stars in 
present day spiral galaxies. This may suggest that 
DLAs are the source of most of material available for 
star formation at high redshift (Lanzetta et al. 1995). 
Later studies have shown that probable optical coun- 
terparts of damped Lya systems at z <^ 1 are galaxies 
of different morphological types (Le Brun et al. 1997). 

Understanding the chemical evolutionary history 
of galaxies is fundamental to the study of galaxy for- 
mation. Pettini et al. (1994, 1997) have performed 
extensive studies on the metallicity of the damped 
Lya systems at high redshift. Their results indicate 
a mean metallicity Zn/H w 1/10 of the solar abun- 
dance with a large dispersion. It was also noted that 
the element Cr is generally less abundant than Zn rel- 
ative to their solar proportion. Since Zn and Cr are 
thought to be made in solar proportions throughout 
the Milky Way based on their abundances in Galactic 
disk and halo stars, the overabundance of Zn to Cr 
found in damped Lya systems was taken as evidence 
that some Cr was depleted from the gas phase by dust 
grains. 

The HIRES W.M.Keck telescope observations have 
opened the opportunity to study Lya absorption lines 
for a large number of elements including N, O, Mg, 
Al, Si, S, Cr, Mn, Fe, Ni and Zn. Lu et al. (1996a), 
Prochaska & Wolfe (1999), Outram et al. (1999) and 
Pettini et al. (2000) have found that the relative abun- 
dance patterns indicate that the bulk of heavy ele- 
ments in these high-redshift galaxies was produced 
by Type II supernovae; there is little evidence for sig- 
nificant contribution from low- to intermediate-mass 
stars or from Type la supernovae. 

However, the current interpretation combining the 
effects of dust depletion and Type II supernovae en- 
richment (Lu et al. 1996a; Prochaska & Wolfe 1999; 
Vladilo 1998) leads to several inconsistencies. In this 
picture the overabundance of Zn relative to Cr or Fe 
found in DLAs is attributed to selective depletion of 
Cr and Fe by dust grains. However, it then turns out 
that such an interpretation is inconsistent with some 
of the other elemental abundance ratios seen in these 



galaxies, in particular N/O, S/Fe, Mn/Fe and Ti/Fe. 
This led Lu et al. (1996a) and Prochaska & Wolfe 
(1999) to conclude that the overabundance of Zn rel- 
ative to Cr in damped Lya galaxies may be intrin- 
sic to their stellar nucleosynthesis. They emphasize 
the need for another physical process to explain the 
damped Lya abundance patterns. 

Furthermore, comparisons of the distribution of 
Fe/H versus redshift (Lu et al. 1996a; Prochaska & 
Wolfe 2000) for the sample of damped Lya galaxies 
with the corresponding relation for the Milky Way 
disk indicate that the damped Lya galaxies are much 
less metal-enriched than the Galactic disk in its past. 
Since there is evidence that depletion of Fe by dust 
grains is relatively unimportant, the difference in the 
enrichment level between the sample of damped Lya 
galaxies and Milky Way disk suggests that DLAs 
are probably not high-redshift spiral disks. A sim- 
ilar conclusion has been reached by Pettini et al. 

(1999) . They argue that the relative abundances of 
heavy elements in DLA systems are consistent with 
a moderate degree of dust depletion that, once ac- 
counted for, leaves no room for the enhancement of 
a elements over iron seen in metal-poor stars in the 
Milky Way. Recently Centurion et al. (2000) have 
found the S/Zn abundance ratio in the relatively high- 
metallicity DLAs to be lower than solar ratio. This is 
contrary to previous assertions that DLAs have been 
enriched solely by Type II supernovae, but it can be 
understood if star formation in these systems pro- 
ceeded at a lower rate than in the early history of our 
Galaxy (Pettini et al. 1999; Centurion et al. 2000). 

It is generally believed that the gas in DLAs is 
mostly in the neutral stage due to the very high opti- 
cal depth beyond the hydrogen ionization limit (e.g., 
Viegas 1994; Prochaska & Wolfe 1996). Therefore, 
it is assumed that all species with ionization poten- 
tials lower than hydrogen are in a singly ionized stage, 
while other species are neutral. With this assumption 
correction factors for ionization are small and the ra- 
tios of column densities of ions are equal to the ratios 
of element abundances. However, Lu et al. (1996a), 
Prochaska & Wolfe (1999) and Wolfe & Prochaska 

(2000) find a good correlation between the velocity 
structure of Al III and singly ionized species. The ion- 
ization potential of Al + is 18.8 eV, i.e. greater than 
that of hydrogen. Therefore, Al +2 is likely present in 
ionized, not in neutral gas. To explain the similarity 
of Al III and other low ionization species line profiles, 
Howk & Sembach (1999) and Izotov & Thuan (1999) 
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proposed that these lines originate in the same ion- 
ized region or in a mix of neutral and ionized clouds, 
and stressed the importance of abundance correction 
for ionization effects. Indirect arguments to consider 
ionized gas in DLAs come from studies of the warm 
ionized medium in the Milky Way and other nearby 
galaxies (Howk & Sembach 1999; Howk, Savage & 
Fabian 1999; Sembach ct al. 2000; Jenkins et al. 
2000). 

Higher ionization species of C IV and Si IV are also 
observed in the spectra of DLAs. However, in the ma- 
jority of cases, the C IV and Si IV profiles differ from 
those of lower ionization species; few cases exist where 
some components of the C IV and Si IV absorption 
lines coincide with e.g. Si II, Al II and Al III (see 
e.g., Lu et al. 1996a; Prochaska & Wolfe 1999; Wolfe 
& Prochaska 2000). The high ionization species are 
thus apparently formed in regions with different con- 
ditions and will not be considered here. In rare cases 
absorption lines from neutral species with ionization 
potentials below that of hydrogen, such as C I and 
Mg I, are observed in DLAs. The existence of such 
lines is a priori not in contradiction with the presence 
of Al III, since neutral species can also be present in 
small quantities in ionized gas (cf. Sect. 3). 

Motivated by the close similarity between the ve- 
locity structures of Al III and singly ionized species in 
DLAs, which suggests the presence of ionized gas in 
the regions where these absorption lines are formed, 
we here consider a simple model that takes these ef- 
fects to first order into account. We examine a model 
of DLAs in which absorption lines of neutral and low- 
ionization heavy element species originate in ionized 
rather than in neutral gas, and we explore the possi- 
ble implications on abundance determinations in such 
systems. As indicated by the presence of high ioniza- 
tion species and the complex structure of the absorp- 
tion line profiles seen in the high resolution spectra, 
the conditions of the absorbing gas must in reality 
be highly variable and of complex nature (see e.g. 
Lauroesch et al. 1996). In the future, both more ad- 
vanced analysis of the detailed line profiles and sub- 
components as well as more sophisticated models al- 
lowing for complex geometries, various ISM phases 
etc. will be necessary to obtain a more accurate de- 
scription of the absorbing systems. It is the hope that 
the present simplified model reflects the main trends 
due to ionization effects. 

The model assumptions are described in Section 2. 
The results from the model calculations and the impli- 



cations on heavy element abundance ratios in DLAs 
are discussed in Section 3. In Section 4 we briefly con- 
sider the application of our model to other objects and 
summarise the main results from this work. 

2. MODELS 

The above discussion shows that the properties of 
the damped Lya systems are yet poorly understood. 
While a significant amount of neutral gas is observed 
in these systems, it is not clear how the regions of ion- 
ized gas are related to those of neutral gas. What is 
the main source of ionization in DLAs - internal stel- 
lar radiation or harder background intergalactic radi- 
ation? Where are absorption lines of heavy elements 
formed - in ionized or neutral gas? Is the metallicity 
of ionized and neutral gas similar or not? Is heavy 
element depletion by dust grains important? 

In the present paper we explore the following 
"multi-component" picture for DLA systems: 1) an 
ionization bounded region illuminated by an internal 
radiation field complemented by 2) a neutral region 
with a lower metal content. These working hypothe- 
ses are adopted to study the main effects intrinsic 
stellar ionizing radiation may have on observed heavy 
element patterns in such conditions. 

To explore the consequences of such a models we 
use the CLOUDY ionization equilibrium code (ver- 
sion c90.05; Ferland 1996; Ferland et al. 1998) to cal- 
culate the ionization state of region 1. For simplicity, 
we assume internal stellar radiation and consider a 
plane-parallel geometry. The outer boundary of the 
slab is determined by the condition that a low electron 
temperature of 2000 K is reached. For the cases dis- 
cussed here (low metallicity and low dust content) this 
boundary condition is a good approximation for the 
outer boundary of ionization-bounded H II regions. 
For the ionizing spectrum, we adopt the spectrum of 
a typical O star. We use an ATLAS line-blanketed 
model atmosphere (Kurucz 1991) with effective tem- 
peratures T cS = 33000 K, 40000 K and 45000K and 
log g = 4.0. The lower temperature is likely more ap- 
propriate for the conditions in the interstellar medium 
of metal-rich spiral galaxies and the local ISM in our 
Galaxy (cf. Sembach et al. 2000). In low-metallicity 
environments the temperature of the ionizing radi- 
ation tend to be higher as evidenced from studies of 
nearby extragalactic giant H II regions (e.g. Stasinska 
& Lcithcrcr 1996). The ionization parameter which 
is proportional to the ratio of the number density of 
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the ionizing photons to that of the electrons is var- 
ied in the range U — 10~ 5 - 10 -2 . We consider the 
metal abundance of 0.1 solar, with relative element 
abundances equivalent to those observed in the low- 
metallicity BCDs (e.g., Izotov & Thuan 1999). In 
particular, C/O and N/O abundance ratios are re- 
duced by factors of 2 and 5 respectively as compared 
to solar abundance ratios. We also adopt for 4 He the 
mass fraction Y — 0.252 (or log He/H = -1.07), which 
is higher than the primordial mass fraction Y p — 0.245 
(Izotov et al. 1999), but is typical for the BCDs with 
the heavy element abundance 1/10 solar (Izotov & 
Thuan 1998). Dust grains with a dust-to-gas ratio 
of 1/10 that in Milky Way are taken into account 
in heating and cooling processes, although the differ- 
ences in ionic fractions derived in the models with and 
without dust grains are small. 

A similar model has recently been considered by 
Howk & Sembach (1999). The major difference with 
their approach is that they consider ionized density- 
bounded regions surrounding the neutral, H I-bearing 
clouds. The boundary of the ionized region is defined 
where the local fraction of neutral hydrogen climbs 
above 10%. With these assumptions the fraction of 
neutral species like O , N° or Ar° is significantly re- 
duced because they always reside near the edges of 
H II region. Additionally, the variation of the local 
ionization fraction of neutral hydrogen at the edge of 
ionized density-bounded region results in large vari- 
ations of some low-ionization species, in particular, 
of Cr + and Fe + . Further resulting differences be- 
tween the model by Howk & Sembach (1999) and our 
ionization-bounded model will be discussed below. 

An important feature of both models is that the 
heavy element abundances in the ionized and neutral 
regions are a priori unrelated. This implies in par- 
ticular that relative abundance determinations with 
respect to hydrogen, as derived from the observed col- 
umn densities, may not be meaningful or are at least 
uncertain. 

3. RESULTS 

3.1. Ionization fractions of elements 

In Table 1 we show the ionization potentials and 
the ionization fractions xf^ of the most important ions 



which are observed or can potentially be observed in 
DLAs, as a function of the ionization parameter U for 
two temperatures T cS = 33000 K and 40000 K of ion- 
izing radiation and the heavy element mass fraction Z 
= O.1Z . On the first line we also show iV(H°+H+), 
the total hydrogen column density across the ionized 
slab. Over the entire range of U, except for the high- 
est value shown here (U = 10~ 2 ), iV(H°+H + ) remains 
significantly smaller than the column density of neu- 
tral hydrogen observed in damped Lya system^. For 
an ionization parameter U ~ 10 -5 hydrogen is essen- 
tially neutral, while at U J> 10~ 3 it is almostly ionized 
as it is evidenced by small neutral hydrogen fraction 
x(H°) (Table 1). Calculations for lower heavy element 
abundances show that the ionization fractions are not 
strongly changed. Therefore, the data from Table 1 
can be used for objects with metallicity [Fe/H] <; -1. 

The distribution of the ionic fractions x(X l ) of 
selected ions as a function of ionization parameter 
U is shown in Figure 1. While the ionization frac- 
tions of all neutral and singly ionized species decrease 
with increasing ionization parameter U, the ioniza- 
tion fraction of Al +2 is increased. Therefore, the ra- 
tio x(Al +2 )/x{A\+) = N(Al III) / 2V(A1 II) is a very 
sensitive indicator of ionization conditions and can 
be used to determine U. Unfortunately, although line 
profiles of Al II and Al III in damped Lya systems 
are essentially similar (Lu et al. 1996a; Prochaska & 
Wolfe 1999), in the many cases the Al II lines are sat- 
urated or marginally saturated. This circumstance 
complicates direct estimates of the ionization param- 
eter U. However, very often, unsaturated lines of Si 
II, Fe II and Al III are observed in the same spectra 
and they can be used to estimate U. For this, we 
assume that the observed heavy element abundance 
patterns at low metallicities are similar to that ob- 
served in low- metallicity halo stars, i.e. [Si/Fe] ~ 0.3, 
[Al/Fe] - -0.4 (Timmes et al. 1995; Samland 1998), 
and adopt that log (Al/Si) Q » log (Al/Fe) Q « -1.0 
(Anders & Grevesse 1989). Then, the expected ratio 
of column densities is equal to 



log 



A(A1 III) 
N(Si II) 



-1.7 + log 



c(Al+ 2 ) 
ar(Si+) : 



(1) 



1 By definition the ionization fraction x(X+ l ) is the ratio of the 
radially integrated number density of ion X+ l to the radially 
integrated total number density of the element X in all ioniza- 



tion stages, i.e. x(X +1 ) = J n(X +l )dr / J n(X)dr, the relative 
column density of the ion X +l . 
2 For T e g = 40000 K and the "typical" value of log U ~ -3.25 
the column density is ~ 40 % of the lower limit which defines 
damped systems. The baryonic mass of low column density 
DLAs could thus be underestimated by up to this amount. 
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log 



AT(A1 III) 
N(Fe II) 



-1.4 + log 



x{Al+ 2 ) 
x(Fc+) ' 



(2) 



The typical values of log [N(Al III) /AT (Si II)] and 
log [N(A\ III)/AT(Fe II)] in damped Lya systems 
are ~ -2.0 to -2.5 and ~ -1.5 to -2.0, respectively, 
which translates to log [x(Al +2 )/a;(Si + )] ~ -0.3 to - 

0.8, [x(Al+ 2 )/x(Fe+)] - -0.1 to -0.6, or log f7 3 

(Fig. 1). This estimated value of ionization param- 
eter is dependent on the uncertainties in abundance 
determination and those for the adopted photoionized 
model and possible depiction of elements into grains. 
In principle, other single ionized species, e.g. Zn II, 
can be used instead of Si II or Fe II. However, Si II 
and Fe II absorption features are more often seen and 
measured in the DLAs, and the atomic data for these 
ions are known with better precision (Howk & Sem- 
bach 1999). The above estimate of U is based on the 
assumption that heavy elements are not present in 
the neutral gas surrounding the region of ionized gas. 
If the neutral gas is significantly enriched in heavy 
elements then an additional contribution to the col- 
umn densities of Si II and Fe II comes from neutral 
gas and the ionization parameter should be increased 
to account for observed N(Al III) /AT (Si II) and N(Al 
III)/iV(Fe II) ratios. 

Before we proceed to discuss the implications of 
this model when applied to observations of damped 
Lya systems it is useful to recall that some neutral 
species with ionization potentials below that of hy- 
drogen can be present in ionized gas as well. For in- 
stance Lu et al. (1996a) have measured logiV(Mg II) 
> 14.38 and logA^Mg I) = 12.44±0.04 in the z damp 
= 0.8598 system toward Q0454+0356 which are con- 
sistent with the predictions of our simplified models 
(Table 1). Their measurements for another Zdamp = 
1.2667 system toward Q0449-1326 provide only lower 
limits logiV(Mg II) > 14.1 and logAT(Mg I) > 12.62, 
also not in contradiction with the ionization model. 
Prochaska & Wolfe (1999) have derived log N(C H) > 
14.6 and logA(C I) < 12.8 in the z damp = 2.140 to- 
ward Q0149+33 which are consistent with predictions 
of ionization models. Hence, the presence of species 
like C I or Mg I does not necessarily mean that they 
reside in neutral gas. 

A useful validity check of our model is the compar- 
ison of the predictions of photoionization models with 
the observed values of [O I/Si II]. In the case of the 
neutral gas model it is expected that [O I/Si II] ~ 0, 
while in the case of our ionization models [O I/Si II] 
= -0.5 - -0.8 depending on T e g and for the range log 



U = -3.5 — 3.0 (Table 1). In the majority of cases 
observations of DLA systems give only lower limits of 

I abundance because of saturation of the O I A 1302 
absorption line. The available data in the literature 
shows that, in general, the lower limits of [O I/Si II] 
are consistent with values expected from the ioniza- 
tion models. However, in some DLAs the lower limits 
of [O I/Si II] are larger than those expected from the 
ionization models. Lu ct al. (1996b) find [O I/Si II] 
> -0.27 for the z = 4.38 absorber towards 1202-0275, 
while Prochaska & Wolfe (1999) find [O I/Si II] = - 
0.17 in the absorber towards Q0951-04. These values 
seems difficult to reconcile with ionization models un- 
less a low ionization parameter is appropriate. How- 
ever, no measurements of Al III have been done for 
these absorbers which would allow to constrain the 
ionization parameter U. 

Recently Molaro et al. (2000) have derived a high 
value of [O I/Si II] = 0.0 - 0.1 from the unsaturated O 

1 and Si II lines in the z = 3.39 absorber towards the 
QSO 0000-2620. For this absorber N(A\ III) / N(Si 
II) = -2.4 is rather high, suggesting an ionization pa- 
rameter \ogU ~ -3.0 to -3.5. However, inspection of 
Fig. 1 in Prochaska & Wolfe (1999) indicates differ- 
ent profiles of Si II A1808 and Al III A1854 in this 
absorber, at variance with other absorbers rendering 
the estimate of the ionization parameter U rather un- 
certain. 

The above discussion suggests that large ionized re- 
gions and abundance discontinuities might be present 
in some DLAs, while other DLAs might be essentially 
neutral. It is likely that DLAs may represent a mixed 
bag of systems with varying ionization conditions. 
While the ionization model does not pretend to ex- 
plain abundance pattern in all damped Lya systems, 
we conclude that in many cases the derived lower lim- 
its of [O I/Si II] do not rule out this model. However, 
we point out that better determinations of [O I/H I] 
are needed to test the importance of the ionization 
effects in the damped Lya systems. 

3.2. Age-metallicity relation and the large 
spread in [Fe/H] 

A very extensive literature exists on the interpre- 
tation of the metallicity of DLA as a function of red- 
shift (e.g. Pettini et al. 1994; Lu et al. 1996a; Lindner, 
Fritze - von Alvenslcben & Fricke 1999; Prantzos & 
Boissier 2000). Among the main issues discussed are 
the question of apparent evolution or non-evolution 
of metallicity with z, the origin of the scatter, and 
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possible observational biases affecting this relation. 
A detailed comparison with such studies is beyond 
the scope of the present paper. Here, we shall only 
indicate the typical metallicity shift expected in the 
framework of our model and mention the main impli- 
cation in terms of age-metallicity relations. 

Pettini et al. (1994, 1999), Lu et al. (1996a), 
Prochaska & Wolfe (2000) noted that the chemical 
enrichment processes in damped Lya galaxies appear 
to be quite inhomogeneous. The [Zn/H] and [Fe/H] 
distributions show a scatter by a factor of <~ 20 with 
the mean value of [Zn/H] which is 0.4 larger than that 
of [Fc/H] for the DLAs at z w 2 (Prochaska & Wolfe 
2000). Furthermore, Lu et al. (1996a) and Pettini et 
al. (1999) noted that the degree of chemical enrich- 
ment in the sample of damped Lya galaxies is, on 
average, considerably lower than the Milky Way disk 
at any given time in the past. This lead them to con- 
clude that the chemical enrichment of damped Lya 
galaxies is different from that in the disk of Milky 
Way. 

A possible explanation for the difference in chem- 
ical enrichment is simply that DLAs are in general 
not necessarily spiral galaxies (Le Brun et al. 1997). 
Alternatively, if in the framework of our model the 
heavy element lines are produced only in the ionized 
region (i.e. the neutral gas has a very low metal con- 
tent), the real Fe/H abundance in the ionized gas is 
several times higher than that derived from the iV(Fe 
II) /N(H I) ratio, because of the non negligible upward 
correction for iron (Table 1). The upward correction 
for zinc is lower, and interestingly it nearly compen- 
sates the difference between [Fe/H] and [Zn/H] dis- 
cussed by Prochaska & Wolfe (2000). Then, the age- 
metallicity relation for damped Lya galaxies is more 
consistent with that for thick disk stars of our Galaxy. 
The opposite situation is of course realised if the neu- 
tral gas has the same metallicity as the ionized gas. 
In this case the observed Fe/H ratio is higher than 
real because in ionized gas x(Fe + ) > x(H°) while in 
neutral gas x(Fe + ) = x(H°) = 1. As we will show 
below, it is, however, difficult to adjust all observed 
heavy element abundance ratios in the frame of the 
model with similar metallicities of neutral and ionized 
gas. 

3.3. Ionization vs Dust Depletion 

There is little doubt that some depletion of heavy 
elements by dust is present in DLAs. However, the 
level of depletion in DLAs is generally smaller than 



that in the disk and possibly in the halo of our Galaxy. 
Pettini et al. (2000) adopt the model of DLA consist- 
ing only of neutral gas and propose for DLAs with 
low level of depletion (systems with typically [Zn/Fe] 
<; 0.3) to correct abundances of refractory elements 
assuming constant depletion factors for all these ele- 
ments. This approach obviously leaves the question 
open how to explain observed abundances in systems 
with more depletion. Besides that, the similarity of 
the line profiles of singly ionized ions and Al III and 
the large iV(Al III)/JV(A1 II) ratio is not addressed. 

Prochaska & Wolfe (1999) have suggested that if 
the overabundance of Si/Fe relative to solar is indica- 
tive of dust depletion, then one might expect a cor- 
relation between [Si/Fe] and [Zn/Fe] with the most 
heavily depleted regions showing the largest Si/Fe and 
Zn/Fe ratios. Indeed, they found such correlation and 
concluded that this appearance is consistent with that 
expected for dust depletion (their Fig. 23). 

Alternatively, as shown in Fig. 2, the regression 
between [Si/Fe] and [Zn/Fe] can be explained by 
ionization effects. Dashed, solid and dotted lines 
in Figure connect theoretical points [x(Si + )/x(Fc + ), 
x(Zn + )/x(Fe + )] in models with Z = 0.1 Z Q calcu- 
lated for various U values and for three values T c g 
= 33000 K, 40000 K and 45000 K of the effective 
temperature characterising the hardness of the ion- 
izing radiation. Note that [x(Si + )/x(Fe + )]~ 1 and 
[a;(Zn + )/a;(Fe + )] _1 are the correction factors to con- 
vert respectively Si + /Fe + to Si/Fe and Zn + /Fe + to 
Zn/Fe and they are independent on the intrinsic 
abundances of Si, Fe and Zn as far as these elements 
do not contribute significantly to the thermal balance 
of the ionized gas. However, the correction factors 
for the considered ions can depend on the helium 
and the total heavy element abundances which reg- 
ulate the thermal balance and the ionization struc- 
ture. As mentioned above, the values of the ioniza- 
tion fractions x are essentially metallicity indepen- 
dent for heavy element mass fractions Z ~ 0.1 - 
0.01 Zq. We point out that the quantities commonly 
(e.g. Prochaska & Wolfe 1999) denoted as [Si/Fe] and 
[Zn/Fe] are actually 

[Si/Fe] = log[JV(Si ll)/N(Fc II)] - log (Si/Fe) Q (3) 

and 

[Zn/Fe] = log[JV(Zn II)/JV(Be II)] - log(Zn/Fe) , 

(4) 

respectively. The total abundance ratios are then dc- 
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fined as 

[Si/Fe] cor - [Si/Fc] - log[x(Si+)/x(Fc+)] (5) 

and 

[Zn/Fc] cor = [Zn/Fe] - log[a;(Zn+)/ a; (Fc + )]. (6) 

We find that both the observed range of x(Al +2 ) 
/ x(Si + ) and [Zn/Fe] are consistent with same ion- 
ization parameter for the ionizing radiation with T e g 
= 33000K and 40000K. Figure 2 shows that if one 
applies the correction for ionization effects for the 
radiation effective temperature of 40000 K, assum- 
ing that after correction [Zn/Fe] cor = 0, the resulting 
[Si/Fe] cor value remains on average above zero with 
a small spread of points (see also Fig. 4), and no 
additional correction for dust depletion is required. 
However, models with harder radiation with T e g = 
45000 K do not reproduce the observed trend. If our 
approach is correct then we obtain important con- 
straint on the hardness of the ionizing stellar radi- 
ation in DLAs from Figure 2. Note that with the 
hard intergalactic background radiation field consid- 
ered by Howk & Sembach (1999) the observed [Zn/Fe] 
vs [Si/Fe] relation cannot be reproduced because of 
their adopted boundary conditions. In this case (see 
their Fig. 1) the correction for ionization of [Zn/Fe] is 
positive, while that of [Si/Fe] is negative. Even if an 
ionization bounded model with the hard intergalactic 
background radiation field is considered, the correc- 
tion for ionization of [Zn/Fe], though smaller, remains 
positive. In our case both corrections are negative. 

In summary, it is possible to reproduce the ob- 
served N(A\ III)/iV(Si II) ratios and the trend seen 
in [Si/Fe] vs [Zn/Fe] diagram assuming ionization by 
a moderately hard stellar ionizing radiation field. In 
principle, no correction for dust depletion needs to 
be invoked, although it might be present in small 
amounts and can easily be taken into account in 
the ionization model without significant changing the 
main features of the [Si/Fe] vs [Zn/Fe] distribution. 

Additional support of the ionization model comes 
from Fig. 3 where we show the dependence of [Al 
III/Fe II] vs [Zn II/Fe II]. Models with T cff = 33000K 
(dashed line) and 40000K (solid line) fit quite well the 
observed points from Prochaska & Wolfe (1999) and 
Lu et al. (1996a) at the values of ionization parameter 
U close to that obtained from Fig. 2. 



3.4. Abundance ratios in ionization bounded 
DLAs 

The above conclusion is broadly consistent with 
expectations for other heavy element abundance ra- 
tios. To illustrate this we apply a "typical" ion- 
ization correction factor to the observational data. 
To do so we adopt the characteristic correction fac- 
tors of the model with gas ionized by the radiation 
with T off = 40000 K. We consider the mctallicity Z 
= 0.1 Zq as representative for the sample of DLA 
systems and adopt log U = -3.25. Then, our cal- 
culations yield the following mean ionization correc- 
tion factors A log[X/H] defined by log[X/H] 

corrected — 

log[X/H] obsorvcd + AlogpC/H]: Alog[Al/Fe] = -0.63, 
Alog[Si/Fe] = -0.31, Alog[S/Fe] = -0.32, Alog[Ti/Fe] 
= -0.12, Alog[Cr/Fe] = +0.16, Alog[Mn/Fe] = -0.01, 
Alog[Ni/Fe] = -0.38, Alog[Zn/Fe] = -0.40. Note the 
small corrections for [Ti/Fe] and [Mn/Fe] which re- 
main almost unchanged while the downward correc- 
tion of [Zn/Fe] is high. Interestingly, it is the same as 
the difference between mean [Fe/H] and [Zn/H] values 
for DLAs at z « 2 from Prochaska & Wolfe (2000). 
The case of nitrogen is discussed separately below. 

The above "typical" correction factors have been 
applied to the sample of damped Lya galaxies from 
Lu et al. (1996a), Prochaska & Wolfe (1999) and Pet- 
tini et al. (2000). Only galaxies with definitely derived 
column densities are included and we do not consider 
objects with lower or upper limits. The original (open 
symbols) and corrected (filled symbols) abundance ra- 
tios are shown in Fig. 4. As pointed out earlier, an 
upward correction of [Fe/H] must also be applied be- 
cause of the large column density of neutral hydro- 
gen. However, this correction is uncertain and does 
not change the relative heavy element abundances. 
Therefore, we simply use [Fe/H] values derived by Lu 
et al. (1996a), Prochaska & Wolfe (1999) and Pcttini 
et al. (2000) keeping in mind that correction might be 
up to of the order of +1 dex. For 4 DLAs of the sam- 
ple the [Fe/H] are not available. We adopt [Fe/H] = 
[Cr/H] for those systems. For [Ni/Fe] we have applied 
an additional upward correction which comes from 
the new measurements of oscillator strengths for Ni 
II from UV absorption lines (Fedchak & Lawler 1999); 
older /-values were used by Prochaska & Wolfe (1999) 
and Lu et al. (1996a). Although the correction is dif- 
ferent for different Ni II transitions, its mean value is 
^0.3 dex, which we adopt for the sake of simplicity. 
We now discuss the behaviour of the element ratios 
of Si, Zn, S, Mn, Cr, Ni, Ti, and Al with respect to 
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Fe. 

Silicon, Zinc: The case of Si and Zn has already 
been treated before (Sect. 3.3). In short, it can be 
seen from Fig. 4 that the corrected [Zn/Fe] values 
scatter around zero, while [Si/Fe] remains positive 
suggesting the production of this element in Type II 
supernovae. 

Sulphur: [S /Fe] shows the same SNII like pattern 
as [Si/Fe]. We note that if the gas in DLA galaxies is 
ionized by softer stellar radiation with T e g = 33000 
K, the downward correction of S/Fe abundance ra- 
tio is 0.17 dex larger and in this case [S/Fe] < 0, 
in contradiction with the production of S in Type II 
supernovae. This difference in S correction factors 
also favours harder ionization radiation for the low- 
metallicity DLAs as compared to that in the Milky 
Way. A new analysis of the S/Zn abundance ratio 
in a sample of DLAs was recently presented by Cen- 
turion ct al. (2000). They find [S/Zn] < in DLAs 
with [Zn/H] ;> -1.0. Both S and Zn are not depleted 
into the dust grains and the correction for dust de- 
pletion does not change S/Zn ratio. Centurion et al. 
(2000) thus conclude that such unusual abundance ra- 
tios suggest that the DLA galaxies are objects with 
low star formation rates such as low surface brightness 
galaxies or dwarf galaxies. Based on our model, the 
correction for the ionization results in the increase of 
both [S/Zn] (Tabic 1) and [Zn/H] and in better agree- 
ment with [S/Fe] for halo stars. If instead the model 
with the hard background intergalactic radiation is 
adopted, the correction of [S/Zn] is negative and it 
increases the difference between [S/Zn] in DLAs and 
[S/Fe] in halo stars. 

Manganese, Chromium: [Mn/Fe] remains nearly 
unchanged because of small correction factor for this 
ratio. The Cr/Fe ratio is overabundant by ~ 0.2 dex 
with respect to the solar value, which seems to be in 
contradiction with the stellar data. This may consti- 
tute a difficulty for our model. Uncertainties in os- 
cillator strengths for Cr II transitions are better than 
10% (Bergeson & Lawler (1993). However uncertain- 
ties in the ionization and recombination cross-sections 
should be examined. Adopting a higher effective tem- 
perature decreases the overabundance of Cr relative 
to Fe, as can be seen from Table 1. The Mn/Fe ratio 
remains underabundant supporting the origin of this 
odd-element in Type II supernovae. 

Nickel: Note that [Ni/Fe] calculated with the new 
oscillator strengths before the ionization correction 
(open circles) are already above zero. This poses 



a problem in the explanation of [Ni/Fe] when dust 
depletion is taken into account, since this correction 
shifts the points upward. Although a large negative 
correction of -0.38 is predicted for [Ni/Fe] in the ion- 
ization model, it is almostly compensated by the up- 
ward shift of points introduced by the new oscillator 
strengths by Fedchak & Lawler (1999). The points in 
[Ni/Fe] vs [Fe/H] diagram scatter around the mean 
[Ni/Fe] value of -0.2. 

Titanium: Prochaska & Wolfe (1999) pointed out 
problems in explaining the Ti/Fe abundance ratio in 
the "dust depletion" model, since Ti is expected to be 
more heavily incorporated into dust grains. Hence, 
the correction for depletion shifts observed [Ti/Fe] to 
larger values. In our model, the correction of [Ti/Fe] 
for ionization effects is negative and small. There- 
fore, its application to Prochaska & Wolfe (1999) data 
leaves Ti moderately overabundant relative to iron 
(filled squares in Fig. 4), in agreement with that ob- 
served in galactic halo stars. 

Aluminium: Our model predicts the largest down- 
ward correction factor of -0.63 dex for Al/Fe abun- 
dance ratio. In the Prochaska & Wolfe (1999) sam- 
ple the [Al/Fe] scatter around zero value. This is 
likely at variance with what is observed in lowest- 
metallicity halo stars where [Al/Fe] scatters below 
zero with mean value of ~ -0.4 (Ryan et al. 1996; 
Timmes et al. 1995; Samland 1998). Our correction 
eliminates this difference and suggests that the odd- 
element Al seen in damped Lya galaxies is likely also 
produced by Type II supernovae. If instead the hard 
intergalactic ionization radiation is applied, the cor- 
rection of [Al/Fe] is <~ +0.1 and shifts points upward 
in the direction opposite to the nucleosynthetic pre- 
dictions. 

Overall, we find that in the framework of our sim- 
plified model of photoionized gas and despite the un- 
certainties in the atomic data, the abundances in DLA 
systems reproduce quite well observed halo star abun- 
dance patterns and lead to a consistent picture. Sev- 
eral problems pertaining to neutral DLA models are 
thus eliminated. In the new model, there is no need to 
take depletion of heavy elements into dust grains into 
account, although small depletion of heavy elements 
might be present. 

In the framework of the proposed model the abun- 
dance ratios of neutral and singly ionized species N 
I / O I, Si II/ S II, Mn II / Fe II (Table 1) depend 
only weakly on ionization effects and can be used as 
a good measure of the total abundance ratios. 
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3.5. On the origin of nitrogen 

We now turn to a discussion of the nitrogen abun- 
dance, which is of particular interest since it is one 
of the elements observable both in metal-poor HII 
regions and DLAs, and because its primary and/or 
secondary origin remains debated (see e.g. Henry, Ed- 
munds & Koppen 2000). 

The heavy element abundance patterns seen in 
DLA galaxies show several similarities with those 
found in local blue compact dwarf galaxies. Izotov 
& Thuan (1999), based on the large sample of low- 
metallicity BCDs, have demonstrated that all ele- 
ments of a process are overabundant relative to iron 
and hence they are produced in Type II supernovae. 

In this respect the N/O ratios seem, at first glance, 
to constitute the major difference. The N/O ratios in 
damped Lya systems appear to be significantly lower 
than those measured in low-metallicity BCDs (cf. Fig. 
5). For the most metal deficient of these objects with 
[O/H] <; -1.3 Thuan, Izotov & Lipovetsky (1995) and 
Izotov & Thuan (1999) have derived a mean value of 
log N/O = -1.60 with a very small dispersion. They 
argue that this value constitutes a lower limit to log 
N/O in BCDs, which is set by primary N produc- 
tion in massive stars. Pettini, Lipman & Hunstead 
(1995) have found an upper limit log N/O <J -1.9 
in one damped Lya system, 0.3 dex below the value 
of Izotov & Thuan. However, within the errors the 
log N/O value derived by Pettini et al. (1995) is still 
consistent with the mean value for BCDs. Further- 
more, Lu, Sargent & Barlow (1998) have compiled 
N/Si abundance ratios in 15 damped Lya galaxies, 
and have set upper limits for log N/O ranging from 
from -1.2 to -2.1. Whereas all these measurements 
constitute upper limits, an actual value of [N/Si] = 
-1.70 has been derived in one system, 1946+7658, by 
Lu et al. (1998). From these DLA measurements of 
N/O and N/Si, Pettini et al. (1995) conclude that 
both primary and secondary nitrogen production are 
important over the whole range of metallicity mea- 
sured in damped Lya systems, while Lu et al. (1998) 
favor the time-delayed primary nitrogen production 
by intermediate-mass stars. 

Can these different scenarios be reconciled ? How 
well are these differences established, especially in 
view of the possible ionization corrections discussed 
above for other elements detected in DLA systems 
? Indeed if our model holds, ionization effects are 
also important for the N/Si abundance ratio. The 



inspection of Table 1 and Figure 1 shows that cor- 
rection factors for the observed A^(N°)/iV(Si + ) ratio 
increase with increasing U. For the typical value log 
U = -3.25 adopted for the above analysis of other 
heavy element abundance ratios (cf. Fig. 4) we find 
Alog[a:(N )/x(Si + )] = 0.72f] 

The observations of N/O in BCDs (Izotov & Thuan 
1999, filled circles) and the DLA data from Lu et al. 
(1998) and Outram et al. (1999) (open circles) assum- 
ing neutral systems are plotted in Figure 5. Over- 
plotted for illustration is the contribution of primary 
and secondary nitrogen production as well as the com- 
bined production by both mechanisms. The "typical" 
ionization correction for N/Si and Si/H^] is shown by 
the arrow in the lower left corner of Figure 5. Tak- 
ing into account such a correction, the vast major- 
ity of measurements for DLAs are found at or above 
the N/O value from the BCDs. This suggests that 
a significant part of the difference in [N/Si] can be 
explained by ionization effects. The typical correc- 
tion brings even the lowest value of [N/Si] ~ -1.7, 
measured in the Zdamp = 2.8448 damped Lya system 
toward Q1946 + 7658, to [N/Si] ~ -0.98 after correc- 
tion for ionization, fairly close to the [N/O] values for 
the BCDs. 

The open square in Fig. 5 shows the N/O abun- 
dance ratio derived recently by Molaro et al. (2000) 
in a DLA with z a b s = 3.3901 from the unsaturated 
O I and N I absorption lines. Potentially, the N/O 
abundance ratio derived directly from N I and O I 
(instead of using Si as a proxy) is the least model- 
dependent. Also it has the advantage of being very 
little dependent on the correction for ionization (see 
Table 1). Most interestingly, the N/O value of Mo- 
laro et al. (2000) is in very good agreement with N/O 
derived in the low-metallicity BCDs. 

In any case, if the typical ionization corrections 
derived from our model apply, we have to conclude 
that the bulk of the present data does not show a sig- 
nificant difference in N/O abundance ratios between 
low-metallicity BCDs and DLA systems implying sim- 
ilar origin of these two elements. We shall now briefly 
summarise recent progress in stellar models regarding 
the nucleosynthetic origin of nitrogen. 

3 In the models of Howk & Sembach (1999) the correction 
Alog[x , (N )/x , (Si + )] is significantly larger (1.48 dex) because 
of the assumption of density-bounded H II region in their case. 

4 Only the correction for Si/H in ionized region is taken into 
account. Larger upward corrections for Si/H should be applied 
if neutral hydrogen region is considered. 
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Interestingly both intermediate-mass and massive 
stars are found to be potential producers of primary 
14 N in the most recent stellar models. Intermedi- 
ate mass stars (Mi n ;ti a i J> 3.5 M Q ) produce primary 
nitrogen when the CN cycle operates in their con- 
vective envelope during the AGB phase (e.g. Ren- 
zini & Voli 1981). This production results from the 
combined effect of the 3rd dredge-up and envelope 
burning, and thus varies both with stellar mass and 
metallicity (e.g. Forestini & Charbonnel (1997) and 
references therein). It is also highly dependent on 
the parameters which intervene in the description of 
AGB stars. Among them, the mixing length and 
the dredge-up parameters which are empirically cal- 
ibrated to reproduce peculiar observed properties of 
different AGB populations (in particular the observed 
luminosity function of carbon stars in the LMC and 
SMC), and which influence the efficiency of the burn- 
ing inside the envelope. In addition, the mass loss 
prescription which remains uncertain (e.g. Vassiliadis 
& Wood 1993) has a crucial impact on the theoretical 
nucleosynthesis by controlling the number of thermal 
pulses undergone by the AGB star and the duration of 
the envelope burning. The most recent models find 
a 14 N production similar to that of Renzini & Voli 
(1981, a = 1.5 case with the largest N production) 
at [O/H] ~ -0.7 but which decreases with increasing 
metallicity (see Marigo 1998; Marigo et al. 1998; Lat- 
tanzio, Forestini & Charbonnel 1999). Intermediate 
mass stars are thus potential primary 14 N produc- 
ers. However, no calculations are yet available in the 
metallicity range of interest for our study. The be- 
haviour at [O/H] < -0.7 remains to be investigated. 

Regarding massive stars it is well known that stan- 
dard stellar evolution models do not produce primary 
nitrogen (e.g. Woosley & Weaver 1995). However, 
considerable advances have been made recently to in- 
clude rotation and its effects on the transport of el- 
ements and angular momentum in the description of 
stellar interiors and evolution (see review of Maeder 
& Meynet 2000). The corresponding new models for 
massive stars, which solve a number of difficulties en- 
countered by "standard" models (cf. Maeder 1995), 
also produce primary nitrogen (Maeder 1999; Langer 
ct al. 1999). In particular the models of Maeder 
(1999) reproduce well the N excess in SMC super- 
giants found by Venn (1999), and suggest that these 
stars could be responsible for primary N production at 
early epochs. Increasing rotational velocities found at 
low O/H (Maeder, Grebel & Mermilliod 1999) could 



favor the primary N production in metal-poor envi- 
ronments. Quantitative yields from these new stellar 
models have still to be awaited. 

Last, but not least, such new yields from interme- 
diate and massive stars will be included in chemical 
evolution models to fully exploit all the observational 
constraints regarding N and other elements. This 
should hopefully solve the question on the origin of 
primary N. Independently of the outcome, the con- 
clusion from our work is that in the light of possible 
ionization corrections for DLAs one cannot presently 
exclude the same primary origin of N in both DLAs 
and nearby metal-poor galaxies. 

4. DISCUSSION AND SUMMARY 

The above analysis shows the importance of consid- 
eration of the ionization effects in the DLAs element 
abundance studies which was pointed out earlier by 
Prochaska & Wolfe (1999), Howk & Sembach (1999) 
and Izotov & Thuan (1999). However, at variance 
with many previous studies assuming that the heavy 
element absorption lines in DLA galaxies are formed 
in the neutral gas, we consider the following simplified 
"multi-component" picture for DLA systems: Region 
1) A plane-parallel ionization bounded region illumi- 
nated by an internal radiation field complemented by 
Region 2), a neutral region with a lower metal con- 
tent. 

Interestingly it appears that the structure of at 
least some DLA systems might be similar to that ob- 
served in nearby low-metallicity blue compact dwarf 
galaxies. For example, the BCD galaxy SBS 0335-052 
showing supcrgiant H II regions over a region of ~ 7 
kpc with ongoing star formation embedded in much 
larger clouds of neutral gas with a column density 
A(H I) as high as 7xl0 21 cm~ 2 , can be considered 
as a local damped Lya system (e.g., Thuan & Izotov 
1997). 

Another interesting parallel is the most metal-poor 
BCD galaxy I Zw 18, which shows a high hydrogen 
column density N(H I) « 3.5 xlO 21 cm" 2 (Kunth et 
al. 1994) thus also qualifying as a local DLA system, 
and for which recent FUSE spectroscopic observa- 
tions were obtained by Vidal-Madjar et al. (2000). 
In contrast with high-rcdshift DLAs the structure of 
I Zw 18 is known better. It consists of two central 
regions of active star formation ionizing the ambient 
gas. Radio observations and the presence of strong 
Lya absorption line indicate that the H II region is 
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embedded in the outer H I cloud. Measurements of 
the column densities of N I, Si II, Ar I and Fc II 
derived from the FUSE spectrum of I Zw 18 were 
kindly made available to us by Alfred Vidal-Madjar 
(2000, private communication). It is striking that the 
column density ratios given in Table 2 (column 2) 
are very different from those observed in the H II re- 
gions (NW and SE components) of I Zw 18 and other 
low-metallicity BCDs also given in the Table 2. We 
suspect that such large differences are due to ioniza- 
tion effects. To analyse the influence of the ioniza- 
tion on the observed column densities we consider a 
spherically-symmetric ionization-bounded H II region 
model with the parameters derived from the observa- 
tions (see e.g. Hunter & Thronson 1995; De Mello et 
al. 1998). In particular, the heavy element mass frac- 
tion in H II region is adopted to be 1/50 Zq, while 
for the neutral gas we assume zero metallicity. 

The results of the calculations are shown in Table 
2, where the observed ratio from FUSE (col. 2) and 
the ionization corrected column density ratio (col. 3) 
are given. Note the good agreement with the abun- 
dance ratios derived from the HST and ground-based 
observations of I Zw 18 and other BCDs, while solar 
ratios do not match corrected ratios. To obtain the 
total hydrogen column density a large amount of gas 
must now be added to that of the ionized gas (cf. Ta- 
ble 2). However, since the heavy element ratios are 
already close to the observed values from the H II 
regions in BCDs, there is little room for additional 
metals in the neutral phase. Hence, this analysis pro- 
vides some support in favor of the basic assumption 
of our model that neutral gas in DLAs systems might 
be more metal-poor as compared to the ionized gas. 

From our simplified model (see Sect. 2) the follow- 
ing main results are obtained (Sect. 3): 

1. Using the observed N(Al III)/iV(Si II) column 
density ratio to constrain the ionization parameter 
U, one obtains almost complete ionization of hydro- 
gen in regions with a column density Y(H°+H+) ~ 
10 20 cm~ 2 (region 1), several times smaller than the 
observed HI column in DLA systems. In the ionized 
region Fe is mainly in the Fe +2 form. Both effects im- 
ply that the true Fe/H abundance in the ionized gas 
must be larger (by up to one order of magnitude) than 
that derived from the iV(Fe II)/AT(H I) ratio. Com- 
pared to [Fe/H] derived from a neutral DLA model, 
such ionized systems have an age-metallicity relation 
in better agreement with that of the thick disk of the 
Milky Way. However, in general, due to the assump- 



tion of a significant metal-deficiency in the neutral re- 
gion (2), abundance ratios with respect to hydrogen 
are a priori not meaningful when the total absorbing 
column is taken into account. 

2. Accounting for ionization effects improves the 
agreement between DLAs and Milky Way stellar thick 
disk measurements. Although not excluded, no signif- 
icant depletion of elements by dust grains is required. 
The majority of heavy element abundance ratios are 
brought into agreement with those expected from the 
production by Type II supernovae. The apparent dif- 
ference between [N/O] derived in local blue compact 
dwarf galaxies and [N/Si] derived in DLA systems 
can be fairly well explained by ionization effects in 
the latter. Thus one cannot currently exclude the 
same primary N origin in both DLAs and metal-poor 
galaxies. 

The above findings lead to a more consistent pic- 
ture for the heavy element pattern in DLAs, and ren- 
der our simple model attractive. 

The main sources of uncertainties potentially af- 
fecting our results arc: poor atomic data (oscillator 
strengths, recombination coefficients etc. for the pho- 
toionization models, affecting e.g. Zn, Cr, and Ni), 
the unknown structure and geometry of the absorbing 
systems, the spatial distribution of ionizing sources 
and their radiation field, and uncertainties related to 
the measurement of column densities. Future investi- 
gations should hopefully be able to improve on these 
issues and to assess the validity of our picture. For 
example, quantitative studies of the line profile struc- 
ture of different ions, and detailed ionization proper- 
ties of individual systems and large samples, should 
be extremely useful. It is the hope that the proposed 
simplified model reflects the main trends due to ion- 
ization effects in DLAs. If correct, our picture of DLA 
systems partly ionized by internal radiation offers a 
clear simplification in the understanding of heavy ele- 
ment abundance ratios in DLAs and their comparison 
with the local Universe. 
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TABLE 1 

Ionic abundance fractions for Z=10~ 1 Z < 



Property IP a r eff =33000K T eff =40000K 







log u — — o.u 


— 4.U 


— O.O 


—o.U 


— z.U 


log U ——o.U 


— 4.U 


— O.O 


— o.U 


— z.U 


iV(H +H+) 




5.24xl0 18 


1.78xl0 18 


4.49 xlO 19 


1.31xl0 20 


1.30xl0 21 


6.38xl0 18 


1.95xl0 19 


4.56xl0 19 


1.29xl0 20 


1.27xl0 21 


log x(H°) 


1.00 


-0.186 


-0.551 


-0.888 


-1.297 


-2.195 


-0.154 


-0.455 


-0.785 


-1.224 


-2.172 


log a; (He ) 


1.81 


-0.008 


-0.029 


-0.040 


-0.047 


-0.044 


-0.097 


-0.321 


-0.581 


-0.897 


-1.193 


log x(C°) 


0.83 


-1.686 


-2.420 


-2.858 


-3.346 


-4.440 


-1.335 


-2.086 


-2.579 


-3.202 


-4.660 


log x(C+) 


0.83,1.79 


-0.009 


-0.004 


-0.006 


-0.016 


-0.118 


-0.024 


-0.037 


-0.111 


-0.279 


-0.795 


log x(C+ 3 ) 


3.52,4.74 


















-5.638 


-4.206 


log x(N°) 


1.07 


-0.205 


-0.595 


-0.947 


-1.376 


-2.305 


-0.184 


-0.506 


-0.852 


-1.317 


-2.333 


log x(N+) 


1.07,2.18 


-0.424 


-0.127 


-0.052 


-0.020 


-0.015 


-0.461 


-0.168 


-0.110 


-0.211 


-0.730 


log x(0°) 


1.00 


-0.176 


-0.534 


-0.870 


-1.278 


-2.183 


-0.145 


-0.441 


-0.770 


-1.210 


-2.182 


log x(Mg°) 


0.56 


-2.290 


-2.548 


-2.974 


-3.577 


-5.013 


-2.242 


-2.449 


-2.711 


-3.218 


-4.580 


log x(Mg + ) 


0.56,1.11 


-0.011 


-0.069 


-0.192 


-0.431 


-1.153 


-0.011 


-0.068 


-0.200 


-0.451 


-1.153 


log x(Al+) 


0.44,1.38 


-0.001 


-0.005 


-0.013 


-0.036 


-0.229 


-0.002 


-0.007 


-0.019 


-0.052 


-0.396 


log x(Al+ 2 ) 


1.38,2.09 


-2.740 


-1.991 


-1.540 


-1.104 


-0.389 


-2.533 


-1.838 


-1.397 


-0.991 


-0.521 


log x(Si+) 


0.60,1.20 


-0.004 


-0.071 


-0.225 


-0.491 


-1.133 


-0.004 


-0.068 


-0.226 


-0.501 


-1.113 


log x(Si+ 3 ) 


2.46,3.32 






-6.061 


-5.351 


-4.147 




-4.570 


-3.423 


-2.612 


-1.489 


log x(P+) 


0.77,1.45 


-0.008 


-0.025 


-0.080 


-0.221 


-0.830 


-0.017 


-0.034 


-0.109 


-0.318 


-0.935 


log x(S+) 


0.76,1.72 


-0.008 


-0.039 


-0.114 


-0.293 


-0.995 


-0.017 


-0.078 


-0.220 


-0.485 


-1.243 


log x(S+ 2 ) 


1.72,2.56 


-1.894 


-1.075 


-0.637 


-0.309 


-0.046 


-1.602 


-0.789 


-0.402 


-0.173 


-0.034 


log x(Ar°) 


1.16 


-0.461 


-0.900 


-1.270 


-1.709 


-2.657 


-0.470 


-0.799 


-1.160 


-1.647 


-2.698 


log x(Ti+) 


0.50,1.00 


-0.022 


-0.157 


-0.372 


-0.698 


-1.516 


-0.023 


-0.153 


-0.380 


-0.738 


-1.619 


log x(Cr+) 


0.50,1.21 


-0.057 


-0.358 


-0.675 


-1.054 


-1.874 


-0.052 


-0.315 


-0.627 


-1.035 


-1.948 


log x(Mn+) 


0.55,1.15 


-0.058 


-0.263 


-0.504 


-0.793 


-1.499 


-0.057 


-0.242 


-0.494 


-0.832 


-1.652 


log x(Fe + ) 


0.58,1.19 


-0.032 


-0.219 


-0.495 


-0.869 


-1.726 


-0.026 


-0.198 


-0.475 


-0.873 


-1.804 


log x(Fe +2 ) 


1.19,2.25 


-1.160 


-0.402 


-0.167 


-0.063 


-0.016 


-1.263 


-0.439 


-0.186 


-0.127 


-0.602 


log x(Co+) 


0.58,1.26 


-0.006 


-0.023 


-0.143 


-0.477 


-1.394 


-0.006 


-0.025 


-0.159 


-0.505 


-1.383 


log x(Ni+) 


0.56,1.34 


-0.006 


-0.016 


-0.113 


-0.427 


-1.350 


-0.006 


-0.019 


-0.138 


-0.474 


-1.346 


log x(Cu+) 


0.57,1.49 


-0.015 


-0.149 


-0.386 


-0.744 


-1.621 


-0.020 


-0.177 


-0.437 


-0.813 


-1.707 


log x(Zn + ) 


0.69,1.32 


-0.006 


-0.005 


-0.012 


-0.034 


-0.234 


-0.017 


-0.063 


-0.170 


-0.375 


-0.888 



ionization potential of the X* 1 ion (first value) and X* ion (second value) expressed in Ryd. 



Table 2 

Abundance ratios in I Zw 18 from FUSE, HST and ground-based observations 



FUSE & Model I Zw 18 (NW) b I Zw 18 (SE) b BCDs c Sun d 



N(H), cm- 2 2.1xl0 21c 7.7xl0 20f 

D, Mpc g • • • 10 

logQ, s~ lh ••• 51.6 

log R in , cm' • • • 20.08 

T off , K j • • • 50000 

N, cm- 3k • • • 100 

f ... o.l 

logN/Si m -0.77 -0.14 ••• -0.09±0.23 ~0.10±0.06 +0.32 

logAr/Si m -1.32 -0.57 -0.61±0.22 -0.74±0.22 -0.72±0.09 -0.95 

logFe/Si m -0.77 -0.42 ••• -0.32±0.24 -0.24±0.12 -0.04 



a Vidal-Madjar et al. (2000), Vidal-Madjar (2000, private communication). 

b Si and Fe abundances are from Izotov & Thuan (1999), abundances of other elements are from Izotov et 
al. (1999). 

c Mean abundance ratios for the BCDs with oxygen abundance 12 + log O/H < 7.60 (Izotov & Thuan 
1999). 

d Anders & Grevesse (1989). 

"Neutral hydrogen column number density. 

'ionized hydrogen column number density. 

s The distance to I Zw 18. 

h Total number of the ionizing photons. 

'Inner radius of the ionized shell. 

j Effective temperature of the ionizing radiation. 

k Number density in H II region. 

'Volume filling factor. 

m Second column: column density ratios TV (NI)/7V (Sill), iV(ArI)/7V (Sill), iV (FeII)/7V (Sill) (dex) from the 
FUSE observations. Other columns: total abundance ratios (dex). 
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log U 

Fig. 1. — Ionic fractions of elements along the line of sight x(X l ) as a function of ionization parameter U for the 
closed plane-parallel H II region models ionized by stellar radiation with (a) T e g = 33000 K and (b) 40000 K. The 
curves for A1+ 2 and Si + are shown by dashed lines. 
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_g 5 I i I i I i I i I i I i I i I I 

0.3 0.4 0.6 0.B 1 1.2 1.4 

[Zn/Fe] 

Fig. 2. — The [Si/Fe] vs [Zn/Fe] relation. The filled circles are observational points from Prochaska & Wolfe 
(1999). Dashed, solid and dotted lines connect theoretical points [x(Si + )/x(Fe+), x(Zn + )/a;(Fc + )] in models with 
Z = 0.1 Zq for varying values of the ionization parameter U (in steps of 0.25 dex) marked by numbers in log scale 
respectively, and for three values of the effective temperature T e g = 33000 K, 40000 K and 45000 K. 
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Fig. 3. — The [Al III/Fe II] vs [Zn II/Fe II] relation. The filled circles are observational points from Prochaska 
& Wolfe (1999) and Lu et al. (1996a). Dashed and solid lines connect theoretical points [x(Al +2 )/x(Fe + )], 
[a;(Zn + )/a;(Fe + )] in models with Z = 0.1 Zq for varying values of the ionization parameter U marked by numbers 
in log scale respectively, and for two values of the effective temperature T e g = 33000 K and 40000 K. 
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[Fe/H] 



[Fe/H] 



Fig. 4. — Heavy element abundance ratios for DLAs from Lu et al. (1996), Prochaska & Wolfe (1999) and Pettini 
et al. (1999). Open symbols are observed abundance ratios while filled symbols arc abundance ratios corrected for 
ionization effects. Squares are data for [Ti/Fe] and [Al/Fe] while circles are for other heavy element abundance 
ratios. The ionization-bounded plane-parallel model with ionization parameter log U = -3.25 and stellar radiation 
temperature T c ff — 40000 K is adopted. 
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-1.5 -1 
[0/H], [Si/H] 



-0.5 



Fig. 5.— [N/O], [N/Si] vs [O/H], [Si/H] diagram. Open circles are [N/Si] vs [Si/H] for DLAs from Lu et al. (1998) 
and Outram ct al. (1999), open square are [N/O] vs [O/H] for DLA from Molaro ct al. (2000), filled circles are [N/O] 
vs [O/H] for blue compact dwarf galaxies from Izotov & Thuan (1999). Dotted lines are diagrams for primary and 
secondary productions of nitrogen, while solid line is dependence [N/O] vs [O/H] when both primary and secondary 
productions of nitrogen are taken into account. Solid line with arrow in the lower left corner shows the value and 
direction of abundance ratio correction for the ionization-bounded plane-parallel model with ionization parameter 
log U — -3.25 and stellar radiation temperature T e g = 40000 K. 
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